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I.  INTRODUCTION 


The  Impact  acceleration  forces  imposed  upon  aircraft  occupants  during 
a  naval  aircraft  ditching  accident  often  result  in  serious  injuries  or 
fatalities,  even  if  there  is  no  direct  impact  with  an  object.  It  is 
important,  therefore,  to  be  able  to  accurately  estimate  the  dynamic  and 
physical  parameters  that  relate  these  forces  to  the  likelihood  of  injury. 

The  Naval  Biodynamics  Laboratory  is  conducting  an  impact  acceleration 
research  program  which  is  focused  on  the  dynamic  response  of  the  human  and 
simian  head/neck  system  as  a  function  of  motion  and  anthropometric  parameters. 

The  accurate  estimation  of  an  injury  likelihood  prediction  model  from 
important  qualitative  and  quantitative  information  should  allow  reliable 
inferences  to  be  made  concerning  injury  probability  as  a  function  of  dynamic 
and  physical  variables.  Previous  technical  reports  [6,  8]  have  discussed 
Injury  likelihood  model  estimation  procedures  based  on  experimental  data. 
Another  technical  report  [4]  described  a  procedure  for  incorporating  a  priori 
information  with  empirical  data. 

This  technical  report  discusses  the  broader  topic  of  incorporating 
various  types  of  auxiliary  information  with  the  dichotomous  "injury" 
observations  to  more  accurately  develop  an  Impact  acceleration  injury 
prediction  model.  The  report  describes  a  general  approach  to  the  beneficial 
use  of  auxiliary  information  in  developing  an  injury  likelihood  prediction 
model  via  probit  analysis.  Continuous  preinjury  auxiliary  information  is 
extremely  important  because  most  of  it  can  be  obtained  without  injuring 
experimental  subjects.  It  is  hoped  that  this  general  approach  will  enable 
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accurate  estimation  without  requiring  an  unethically  large  number  of 
dichotomous  injury  observations. 
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II.  BACKGROUND 


The  auxiliary  information  considered  in  this  technical  report  can  be 
classified  into  two  main  types: 

1.  prior  information 

2.  auxiliary  empirical  information. 

The  prior  information  can  consist  of 

(i)  a  priori  estimates  of  some  (possibly  vector-valued)  function  of 
the  model  parameters 

(il)  a  priori  knowledge  in  the  form  of  model  parameter  equality  con¬ 
straints. 

In  the  Injury  probability  prediction  model,  auxiliary  empirical  information 
could  be  obtained  in  the  form  of  an  experimental  "side  effect",  e.g.,  evoked 
potential  response  data. 

The  Injury  likelihood  model  considered  in  this  report  is  to  be  developed 
under  the  assumption  that  the  injury  tolerance  is  normally  distributed.  This 
assumption  seems  to  hold  true  in  many  probability  prediction  settings. 
Furthermore  the  normal  distribution  provides  the  most  tractable,  fundamental 
distribution  model  for  incorporating  auxiliary  empirical  information.  The 
analysis  of  dichotomous  response  prediction  problems  with  normal  tolerances 
is  referred  to  as  "probit"  analysis.  In  previous  technical  reports  [3,  4,  5, 
6,  7,  8],  the  tolerance  was  assumed  to  have  a  logistic  distribution.  This 
assumption  was  made  for  purposes  of  mathematical  tractability  since,  in  those 
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reports,  the  assumption  of  normal  versus  logistic  tolerances  made  little 
difference  in  injury  probability  prediction. 

The  basic  multipredictor  probit  model  is  based  on  the  concept  of  a 
subject  having  a  "tolerance”  for  a  stimulus  or  set  of  stimuli,  e.g.  peak 
sled  acceleration.  The  vector  of  predictor  variables  is  represented  by 

X*  -  (x1,  x2,...,  x^, 

where  at  least  one  of  the  elements  of  x.  represents  a  particular  level  of 
stimulus.  Some  of  the  elements  of  x  may  also  represent  attributes  of  the 
subject  exposed  to  the  stimuli. 

The  problt  model  is  derived  from  the  notion  that  the  stimuli  result 
in  some  nonobservable  real  or  hypothetical  quantitative  response,  t.  The 
observable  dichotomous  response,  y,  (e.g.,  injury  or  no  injury)  is  deemed 
to  occur  if  t  falls  below  some  critical  value,  which  without  loss  of 
generality  may  be  taken  as  zero.  If  the  mean  of  t  is  a  linear  function  of 
jc,  then  we  may  write 

t  -  x'B  +  T 
- T 

where  T  is  a  normal  random  variable  with  zero  mean  and  standard  deviation  a ■  , 
Here  t  is  the  tolerance  of  the  subject,  a  normally  distributed  random  vari¬ 
able  with  mean  x’j^.  and  a  standard  deviation  equal  to  0T>  For  a  given  set 
of  stimuli  the  dichotomous  variable  y  is  set  equal  to  1  if  a  response  occurs 
and  is  set  equal  to  0  otherwise.  Thus, 

Pr(y  -  l|x)  ■  Pr(t  <  Ojx)  -  ♦(-x.'jS j/aT) 

where  $  is  the  standard  normal  cumulative  distribution  function. 

With  the  reparameterization 
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the  following  expression  results: 

Pr(y  -  l]x)  ■  (2) 

The  equation  given  by  (2)  represents  the  basic  multipredictor  probit  model. 

For  reference  purposes,  the  random  variable  y  will  be  denoted  by  y^  for  the 
remainder  of  this  report.  Model  (2)  can  also  be  written  in  terms  of  the 
actual  dichotomous  observations  as: 

yx  -  $(x'B1)  +  ex  (3) 

where  E^)  -  0,  VarCe^  -  p(l  -  p),  p  - 

In  this  technical  report  it  is  assumed  that  for  each  vector  of  stimuli, 
a  corresponding  continuous  auxiliary  response  can  be  measured.  This  response 
could  be  interpreted  as  a  side  effect  due  to  the  stimuli.  It  is  also  assumed 
that  the  relationship  between  the  stimuli  and  the  side  effect  can  be  explaineed 
by  a  linear  regression  model  of  the  form 

y2  ■  x'i2  +  e2,  (4) 

where  y^  is  the  observed  side  effect.  Further  it  is  assumed  that  x  and  £2 
have  a  bivariate  normal  distribution  with  correlation  coefficient  p  . 

It  should  be  noted  that  in  the  continuous  pre injury /dichotomous  injury 
setting  that  p  is  negative.  This  is  so  because  the  subject's  tolerance 
should  be  inversely  related  to  the  level  of  side  effect.  In  other  words, 
as  the  levels  of  stimuli  increase,  the  side  effect  Increases  and  the  subject's 
tolerance  for  the  increase  levels  of  stimuli  decreases.  The  experimental 
data  is  then  modeled  by: 
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i  -  f  (6)  +  e. 


where 
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The  error  structure  has  the  following  properties: 


(a) 


'li 


'2i 


are  independent  vectors  for  i  ■  1,. 


(b)  Eli  and  £2^  are  correlated  for  each  i,  i 


(c)  E(e)  -  0 
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(d)  Var(e)  =  V  = 


*1 


V 

“2 


where  V.  ■  Var 


?i  - 


li 


'2i 


P±(l  -  p±)  Ci 


for  i  a  • i  n  } 


-  -po2(2ir)iSexp{-(x^^1)2/2}. 


(Derivation  of  the  expression  for  is  given  in  the  appendix.) 

As  an  aid  to  both  estimation  and  computation,  a  priori  information 
about  the  parameters  in  model  (5)  can  be  quite  useful.  Any  good  prior 
estimates  of  the  model  parameters  or  some  function  of  them  not  only  help 
obtain  accurate  final  model  parameter  estimates  but  also  aid  in  the  con¬ 
vergence  of  the  estimation  procedure  which  must  be  performed  iteratively. 
Further,  prior  information  in  the  form  of  proper  parameter  constraints  can 
significantly  reduce  the  variance  of  the  constrained  parameters.  As 
presented  in  this  report,  the  prior  information  about  the  model  parameters 
need  not  be  complete.  In  other  words,  prior  information  on  each  individual 
parameter  need  not  exist. 

With  a  few  reasonable  assumptions,  prior  estimates  that  are  in  the 
form  of  some  function  of  model  parameters,  can  be  modeled  in  a  fashion  that 
lends  itself  to  computation  by  least  squares.  For  an  example  of  this  in 
the  linear  least  squares  setting,  see  [9],  To  see  how  this  prior  information 
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can  be  modeled  in  a  general  nonlinear  least  squares  setting,  denote  a 
set  of  generic  model  parameters  by  a  q  x  1  vector,  a.  Let  r  be  a  d  x  1 
vector,  (d  q),  that  is  an  unbiased  estimate  of  some  function  of  a,  g(a) , 
say.  Then  one  can  write 


£  *  &(a)  +  v 


(6) 


where  E(v)  ■  0. 

Assume  further  that  Var(v)  °  T,  where  is  known  or  can  at  least  be 
estimated,  and  that  the  prior  information  embodied  in  r_  is  stochastically 
independent  of  Thus  it  follows,  trivially,  that  £  and  £  are  stochastically 
independent  and  that 


Var 


1 

0 


Further,  if  a  maximum  likelihood  estimate  of  a  is  desired,  it  is  convenient 
to  assume  that  _v  is  multivariate  normal  with  zero  mean  and  covariance  matrix 

I- 

One  may  wish  to  impose  restrictions  across  the  elements  of  a  ,  The 
most  elegant  way  to  represent  these  restrictions  is  by  reparameterization. 
Let  B  be  an  s  x  1  vector  of  new  parameters  relating  a  to  B  according  to 


a  **  h(0)  . 

It  is  assumed  that  s  £  q  .  Thus  under  these  restrictions,  the  "prior 
estimate"  model  in  (6)  can  be  written  as 


r  *  &(h(0)  )  +  v 


(7) 


-8- 


and  h(0)  is  also  substituted  for  a  into  the  model  for  the  empirical  data 


Parameter  constraints  can  also  occur  in  the  form  of  inequalities.  For 
example,  suppose  it  is  known  that 

g^a)  £  0 


Sj  fa)  1  0 


for  some  functions  g^,...,  g^ 


Sometimes  these  inequality  constraints 


can  be  converted  to  equality  constraints  from  results  of  nonlinear  program¬ 


ming  theory,  e.g.,  the  Kuhn-Tucker  conditions.  (See  [1].) 

If  any  of  these  inequality  constraints  cannot  be  removed,  then  formal 
techniques  of  nonlinear  programming  need  to  be  introduced  to  estimate  the 
parameters.  In  such  a  case  the  reader  is  referred  to  [1],  General  nonlinear 


programming  procedures  are  not  considered  in  this  technical  report  as  they 
could  not  be  directly  performed  by  the  basic  nonlinear,  least  squares  statis¬ 
tical  packages  available. 


For  Che  general  estimation  problem  considered  in  this  report,  nonlinear 
weighted  least  squares  estimates  are  most  often  easier  to  obtain  than  are 
the  maximum  likelihood  estimates.  This  is  especially  true  when  there  are 
parameter  constraints  linking  3^  and  J3^.  The  likelihood  function,  however, 
offers  an  enlightening  insight  into  the  relationship  between  the  dichotomous 
injury  data  and  the  continuous  preinjury  data. 

In  this  report,  therefore,  a  theoretical  analysis  of  the  empirical 
information  is  presented  via  the  likelihood  function.  A  future  report  will 
discuss  the  nonlinear  weighted  least  squares  approach  from  a  computational 
point  of  view.  No  real  gaps  in  this  approach  to  the  estimation  problem  will 
occur,  at  least  asymptotically,  since  maximum  likelihood  estimates  and 
weighted  least  squares  estimates  are  essentially  the  same  for  moderate  to 
large  sample  sizes  if  the  data  are  sampled  from  a  distribution  belonging  to 
a  regular,  exponential  family  as  is  the  case  here. 

The  estimation  theory  discussed  here  centers  on  just  the  empirical 
information  in  order  to  keep  the  presentation  from  becoming  too  lengthy. 
Afterwards,  the  incorporation  of  the  a  priori  information  can  be  developed 
in  a  straightforward  manner. 

Denote  the  random  sample,  conditional  or  x^,...,  x^,  by 
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(8) 
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where  as  before 


yil’“‘ 
^21*  *  *  * 


5.^1  •  •  •  * 


,  y  are  the  dichotomous  injury  observations 
In 

,  y2n  are  the  continuous  preinjury  observations 
are  the  corresponding  predictor  variables. 


To  derive  the  likelihood  of  (8),  first  consider  the  joint  distribution 
of  each 


i  *  1, . . . ,  n  , 


2i 


where  as  before 


t^ .  tQ  are  the  unobservable  injury  tolerances. 


Each 


is  by  assumption,  bivariate  normal  with  mean 


x'B 

-1% 


'2i 


U ihi 


and  covariance  matrix 


POt02 


paTa 


x  2 


I  Ci 


For  each  i,  let  f  (t  ,  y9  )  be  the  joint  density  function  of 

2  ly2i 

Then  f£  y  (t^,  y^)  can  be  written  as 


£t(tl|y2i>ty2<y21) 


where  1^21^  *9  t*ie  un*var*ate»  normal  density  of  t^  conditional  on  y2^ 
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It  follows  that 


and  fy  (y^^)  Is  the  univariate,  normal  density  of  y^. 
from  *  t*xe  Pro^ability  distribution  ef  y^  conditional  on  y2^ 

is  obtained.  It  is  generally  known  from  normal  distribution  theory  that 
tily2i  **a8  a  ®esn  equal  to 


ij£T  +  pfl(y2i  "  2EJi2) 
a2 

and  a  variance  of 

<\U  -  P2)  . 


Now  Pr(yil  -  1|  y21)  -  Pr (T^  <  0|xif  y21) 

“  +  P^r(y2i  -  x!^2)]/a  [1  -  p2]*) 

a2 

m  pV^f-VV  -  P^y2i  -  2Lii2)/02>]>- 

However,  -  -j^/cr^  from  (1),  so  that 

Pr(yu  -  ljx^  y21)  «$((!-  P2)-*5!^  -  p((y2i  -  x^/c^)]). 


For  brevity's  sake  let 

Pi  -  Pr(yu  =  ll^,  y2i)  . 


Then 


yliM 

pi  (1“Pi) 


a  -  yu) 


is  the  probability  distribution  of  y,.  conditional  on  x.  and  y  ,  so  that 

li  — ' 


2i 


Pi^U  -  Pl)(1  “  yii)fy2(y2i) 
is  the  joint  distribution  function  of 


yli 


,  y2i 


conditional  on  x. . 

—1 


Thus  the  likelihood  of  (8)  is 


(9) 


n 

n 

i-1 


rli 


a  -  p/1 '  yu>£s,2<i'21) 


where  p^  is  as  previously  defined  and 

fy2(y2i)  “  <2ir)“JSff2’lexi,[~,s<(y2i  "  5il2)/02)2]* 

Now  after  some  reparameterization  a  surprising  result  emerges.  The 
maximum  likelihood  estimates  of  an^  a2  not  depend  upon  y^n  • 

They  are  exactly  the  same  estimates  that  would  be  obtained  from  computing 
the  maximum  likelihood  estimates  of  JB2  and  02  based  on  only  the  preinjury 
data.  To  make  this  clear,  consider  the  fact  that  may  be  expressed  as 


p-t  -  (*Ji*  +  yy21) 


(10) 


where  J3*  -  (1-  p2)”**^  +  PjL2/°2) 
and  y  -  -(1  -  p2)~\/o2  . 

Due  to  the  parameters  and  p  in  (10),  g*  and  y  are  not  explicit 

functions  of  ^  ar>d  02  alone.  The  likelihood  function  in  (9)  then,  can  be 

maximized  with  respect  to  j3  and  02  independently  of  any  of  the 

p^  ^(1  -  p  )^  yli^  functions  in  which  only  the  jj*  and  y  parameters 

appear.  Thus,  the  maximum  likelihood  estimates  of  and  do  not  depend 

upon  the  dichotomous  injury  observations.  This  means  that  the  y  , ...»  y 

11  7  In 

observations  cannot  be  used  to  improve  the  parameter  estimates  of  and  02. 
This  is  not  too  serious  however,  as  these  are  simply  auxiliary  parameters. 
The  maximum  likelihood  estimates  of  £*  and  y  though,  depend  upon  both 


the  y 


11’ 


y^  and  the  y2^,...,  observations.  To  see  this,  let  8*,  y, 


j$2  and  02  be  the  maximum  likelihood  estimates.  Then  one  can  use  the  equations 
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i  A  ^ 

in  (10)  to  solve  for  the  maxinun  likelihood  estimates'1'  3.  and  p;  they  are 

A  M  i 

P  -  -Y<J2/(Y  Pj  +  *> 

A  A  A-J  U  A  A  A  _|_ 

and  -  £*(1  -  pV  -  Pj^ 

A  further  result  of  interest  is  that  the  maxinun  likelihood  estimates 
of  8^  and  p  are  functions  of  the  y^, y2n  observations  only  through  the 
standardised  residuals  from  the  (maximisn  likelihood)  linear  regressions  of 
the  y^'s  on  the  x^'a  .  This  is  shown  as  follows.  Let 

A  A 

m  (y2^  "  — /q2»  ^  "  !*•••»  n 

be  the  standardised  residuals.  Mow  consider  again  the  likelihood  function 
expressed  by  equation  (9).  For  notational  brevity  write  this  function  as 

L(8^  i  P »  .8 2»  ®2)  •  (11) 

The  function  in  (11)  can  be  naxlmised  with  respect  to  (j5^,  p,  8.2,  o2)  by 
maximizing  the  following  function  of  8^  and  p  : 

A  A 

max  L(8lt  p,  £2,  o2)  -  L(8lt  p,  6^,  o) 

I2*°2  “ 

However , 

*  n  a  Vli  *  (1  -  y  )A 

l(ax.  p.  i2,  o2)  -  Pl  Al(i  -  Pl)v  11  fy2(y2i)  (12> 

where  Pj  -  *((1  -  02)~l5[x^_^1  -  ps1]) 

*nd  ^y(Y2i)  “  lexpf-<y2i  ‘ 

Z  I 

8,1  and  p  are  the  maximum  likelihood  estimates  by  the  invariance 
princiPle  of  maximum  likelihood  estimation.  (See  [2].) 


A  A 

since  ■i  -  (y21  -  x|J32)/a2 

Thus,  the  max  1mm  likelihood  estimates  of  and  p  can  be  based  upon  the 
observations 


(13) 


through  the  likelihood  function 

n  p  yii(i  -  it)(1  '  yii}  (l 

i-l  1 

So  it  is  now  evident  that  the  maximum  likelihood  estimates  of  and  p 
can  be  based  upon  (13)  through  (14)  or  based  upon 


yln 


t  •  •  •  9 


(15) 


y2l/ 


r2n 


through  the  likelihood  function 


n  p/11U  -  pt) 

i-l  1  1 


(1  -  yu) 


(16) 


where  pi  -  $(x{g»  +  yy2i)  »  1*1 . . 

However,  the  likelihood  function  in  (16)  has  the  disadvantage  that  the 
continuous  preinjury,  y2i,  observations  are  correlated  with  the  predictor 
variables.  This  could  lead  to  a  problem  of  multicollinearity  in  the 
estimation  of  j5*  and  y  . 
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Asymptotically,  the  nonlinear  weighted  least  squares  estimates  of 
6^  and  p  are  equivalent  to  their  maximum  likelihood  estimates.  The 
nonlinear  weighted  least  squares  estimates  are  derived  from  the  model 

yli  “  *tu1<B1.  p>]  +  £±.  1-1,...,  n  (17) 

where  p)  -  (1  -  p2) -  ps1>,  i  -  1,...,  n 

and  E(f^)  -  0,  Var^)  -  p1(l  -  p1),  p±  -  $[Ui(31,  p)J. 

Computational  aspects  of  this  model  will  be  addressed  in  a  future  report. 
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APPENDIX 


This  appendix  contains  a  detailed  derivation  of 


and  Var(e2iyli)  1  -  1,...,  n  . 

First  consider  the  joint  distribution  of  e ^  and  t^.  By  assumption, 
and  t^  have  a  bivariate  normal  distribution.  For  brevity,  let  the  mean 


of  t^  be  u  .  Now  E(e21y1;l)  equals 


00  00 

f  f  E  I  (t  )f  <e2l,  t  )d£  dt 

-oo  ^  211  i  e,t  21  i  21  1 


where  I±(t1>  - 


0  If  t±  >  0 
1  if  t  <  0 


and  f^  t(e2i»  21s  the  Joint  density  of  e2i  and  t^. 

Thus, 

E(e2iyli)  "  _£  _{e:2ife(E2ilti)ft^ti)de2idtl 

where  fg^jJ^)  is  tl>e  norma-*-  density  of  conditional  on  t^ 
and  ffc(t^)  is  the  marginal  density  of  t^. 

By  integration  with  respect  to  e2^,  one  obtains 


E(e2iyli)  “  7  E<e2ilti)ft(ti)dti< 


(1) 


Now,  since 


-  P^t.  .  > 

Ai 


aTvwi 


(2) 
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by  substitution  of  (2)  into  the  right  hand  side  of  (1),  one  obtains 


E^e2iyli^  “  po2  f 


°/ti- 


i j  ft(t1)dt1  . 


Thus,  -Ji  /O 

Ti  T  h 

E^e2iyli>  "  p°2  ^  exPt-(v^/2)]/[(2Tr)  v1]dvi 


where  v 

i 


<*i 


-pa  <K-u  /a  ) 

2  T±  T 

2  *5 

u  )/a  and  <|)(h)  =*  exp(-h  /2)/(2ir)  . 
Ti  1 


However,  =  x^$  and  3^  ”  -jS/a^  . 

Therefore,  E(e  y^)  *  -po^OE^)  •  Because  E(e2le11)  =  E<e2iyli^’ 
this  is  the  value  of  which  appears  in  section  II  of  this  report. 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  of  This  PAGE  r<Fh«n  Oar.  Enttred I 

!~  REPORT  DOCUMENTATION  PAGE 


I  REPORT  MOMBER 

112-6/ 


|Z.  GOVT  ACCESSION  NO. 


READ  INSTRUCTIONS 
BEFORE  COMPI.ETINCi  FORM 
3  RECIPIENT’S  CATALOG  NUMBER 


A.  TITLE  (and  Submit) 

)  Ji  GENERAL  STATISTICAL  APPROACH  FORCING  ^JXILIARl 
,  INFORMATION  IN  THE  DEVELOPMENT  OF  AN  JMPACT 
|  ACCELERATION  INJURY" PREDICTION  110DEL  / 


John  J.  /Peterson-^^t:  Dennis  E. /Smith 


S^JVPE  OF  REPORT  ft  PERIOD  COVERED 

(^Technical  ^epeit  t  ± 

(.  PERFORMING  ORG.REPSWTNuMBER 
S.  CONTRACT  OR  GRANT  NUM8ERf*J 

\  NOO014-79-C-0128/  / 


».  PERFORMING  ORGANIZATION  NAME  ANO  AOORESS 

Desmatics,  Inc./ 

P.  0.  Box  618 


/'  m 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  «  WORK  UNIT  NUMBERS 

NR  207-037 


IT  CONTROLLING  OFFICE  NAME  ANO  AOORESS 

Office  of  Naval  Research 
Arlington,  Va  22217 


X2.  -RE.PORT  DATE 

1  Auguafc-J.980' 


4.  MONITORING  AGENCY  NAME  ft  ADORESSfff  dtilarent  from  Controlling  Office)  |  1$.  SECURITY  CLASS,  (of  thla  report) 


/f  :•  r  •  - 


Unclassified 

ISa.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


I  IS.  DISTRIBUTION  STATEMENT  (ol  ihl l  Report) 


Distribution  of  this  report  Is  unlimited. 


I  7.  DISTRIBUTION  STATEMENT  (of  the  abatract  ^.r  farad  In  Block  20,  If  dl/farant  from  Report) 


If.  SUPPLEMENTARY  NOTES 


1 19.  KEY  WOROS  (Con tin 


ary  and  Identify  by  block  number) 


Impact  Acceleration  Injury  Prediction 
Use  of  Auxiliary  Information 
Statistical  Model 

20  ABSTRACT  (Continue  on  reaerae  aide  It  neceaaary  end  Identity  by  block  oumborj 

This  report  discusses  general  procedures  for  simultaneously  incorporating 
various  sources  of  auxiliary  information  into  an  Impact  acceleration  injury 
prediction  model.  The  sources  of  auxiliary  information  considered  are  sup¬ 
plemental  continuous  empirical  data  and  a  priori  knowledge  in  the  form  of  model 
parameter  estimates  and  constraints. 


EOITION  of  I  NOV  *S  1$  OBSOLETE 


_ UNCLASSIFIED  -  If 

SECURITY  CLASSIFICATION  of  This  PAGE  (•»*»"  Daia  EmnM 


■11- LiVy. 


> 


